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Abstract 

The M^-pair production in association with a massive (anti)bottom jet is not only an 
important background to a number of interesting processes, such as the single top production 
associated with a W boson, but also a potential background to new physics searches. We 
present the calculations of the total and differential cross sections for the W + W~ + b(b) 
jet productions at the LHC up to the QCD next-to-leading order (NLO). Our results by 
adopting the QCD NLO contribution collection scheme-I show that the K factors can be 
1.66 and 1.21 with the inclusive and exclusive two-jet event selection schemes respectively, 
when we set mn = 120 GeV, a. = raw + nib/ 2 and take the constraints of p T b ^ > 25 GeV, 
\yub)\ < 2.5 for b(b) jet. We find that the stabilization of the theoretical prediction for 
the integrated cross section for the pp — > W + W~b(b) + X up to the QCD NLO requires a 
veto on a second isolated hard jet and the inclusion of the QCD NLO contribution from the 
W + W~bb(bb, bb) production with the final two b(b) quarks being merged as one jet. 



PACS: 12.38.Bx, 13.85.Lg, 14.65.Ha, 14.70.Fm 



1 



I. Introduction 



At the LHC, the single top associated production, pp — >• bg — >• Wt+X, offers a unique possibility 
of the direct measurement of the entry V#, of the Cabibbo-Kobayashi-Maskawa quark-mixing 
matrix (CKM), allowing nontrivial tests of the properties of this matrix in the standard model 
(SM) [U El [3] . By this channel we can also study the W — t — b vertex, and test precisely the V-A 
structure of the charged current weak interaction of the top quark by looking at the polarization 
of this quark[31 0]. Furthermore, such a production channel could be interesting in hunting 
for new physics beyond the SM. The new physics may manifest itself via either loop effects, 
or inducing non-SM weak interactions to introduce new single top production channels [5]. The 
produced top quark of the pp — > bg — )• Wt + X process subsequently decays mainly into a W 
boson and a bottom quark, therefore, the observed final states of the single top-quark production 
in association with a W boson are mostly W + W~b where the two W bosons can be reconstructed 
from their leptonic decay products. The similar analysis also applies for the pp — > W + W~b + X 
process. 

Besides, the Higgs boson can be produced in association with a high pt jet. If the SM Higgs 
boson is in the intermediate mass range, the pp -> #°+jet -> W + W~+iet process is expected to 
be a discovery channel for the Higgs boson, particularly if the Higgs-boson mass is very close to 
the threshold for the VF-pair production. B. Mellado, et al, conclude in their paper [B] that the 
Higgs signal significance can be improved by considering also the final states which consist of the 
leptonic decays of the W pair, and at least one additional jet. They suggest that the Higgs boson 
associated with only a single jet at large rapidity is requested as a signal. They find that the 
backgrounds can be reduced substantially and the sensitivity of experiments can be enhanced 
in the search for the Higgs boson via the W + W~ + jet production channel. The QCD next-to- 
leading order (NLO) corrections to the process pp — > W + W~j + X were studied in Refs.[71 [5]. 
In these papers the initial quark q ^ b,b was assumed, and the subprocesses with j = b,b in final 
states are excluded for the reason that such processes can be excluded by anti-6 tagging or should 
be seen as continuations of W~t or W + t production. Since the process pp — > W + W~b(b) + X 
at the LHC is an essential background of both the single top production and the single Higgs- 
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boson production channels, the investigation of it can serve as a complementary work for those 
presented in Refs.[71[S]- Furthermore, the process pp — > W + W~b(b) +X is interesting in its own 
right, since IF-pair production processes enable a direct precise probing of the nonabelian gauge 
boson couplings. To sum up, improving the precision of the predictions for the W + W~b(b) 
production is necessary, since W + W~ + 6(6)-jet production delivers not only the background to 
a number of interesting processes, but also a potential background to new physics searches. 

In this work, we perform the calculations for the QCD NLO corrections to the pp — > 
W + W~b(b) + X process. Within this work, the mass of bottom quark is retained in all the 
partonic processes. The paper is organized as follows. We present the details of the calculation 
strategies in Sec. II. The numerical results and discussions are given in Sec. Ill, and finally a 
short summary is given. 

II. Strategies in calculation 

We neglect the quark mixing between the two light and the third massive generations (i.e., 
V u b = V c b = Vtd = Vts ~ 0), and take the u, d, c, s quark being massless (m u = ma = m c = 
m s = 0) and the bottom-quark being massive. We use the five-flavor scheme (5FS) in the 
leading-order (LO) and QCD NLO calculations. The difference between adopting the five-flavor 
scheme and four-flavor scheme is the ordering of the perturbative series for the production cross 
section. In the four-flavor scheme the perturbative series is ordered strictly by powers of the 
strong coupling a s , while in the 5FS the introduction of the b quark parton distribution function 
(PDF) allows to resum terms of the form a™ln(^ 2 /m^) m at all orders in a s . The calculations 
are carried out in the 't Hooft-Feynman gauge. The FeynArts 3.4 package[9] is adopted to 
generate Feynman diagrams and convert them to the corresponding amplitudes. FormCalc 
5.4 programs |10j are implemented to simplify the amplitudes. As we know the cross section 
for the gb — > W + W~b partonic process in the SM should be the same as that for its charge 
conjugate subprocess gb — > W + W~b, and the luminosity of the bottom-quark in proton is the 
same as that of the antibottom quark. Therefore, the production rates for both the W + W~b 
and the W + W~b productions at the LHC are identical. In the following sections we present 
only the analytical calculations of the related partonic process gb —> W + W~b and the parent 
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Figure 1: The LO Feynman diagrams for the partonic process gb — > W + W b , where the 
internal top-quark in Figs. 1(4, 6) and the Higgs boson in Figs. 1(1, 7) may be on mass shell. 

process pp — > W + W~b + X unless otherwise indicated. 
II.. 1 Leading order cross sections 

We denote the concerned partonic process as g(pi) + 6(^2) — > W + (ps) + W~{p^) + 6(ps). There 
are nine leading order (LO) Feynman diagrams for this partonic process shown in FigHJ There 
Figs. 1(1-4) and Figs. 1(5-9) are the s-channel and t-channel diagrams for the partonic process, 
respectively. The LO cross section for the partonic process gb —> W + W~b is obtained by using 
the following formula: 



where d&3 is the three-body phase-space element, and p\ is the momentum of the initial gluon 
in the center-of-mass system. The integration is performed over the three-body phase space of 
the final particles W + W~b. The summation is taken over the spins and colors of the initial and 
final states, and the bar over the summation indicates averaging over the intrinsic degrees of 
freedom of initial partons. 

In the LO calculations, the internal top quark and Higgs boson are potentially resonant. 
We introduce the complex mass scheme(CMS) [TT] to deal with the internal resonant top quark 




(2.1) 
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and Higgs boson. In the CMS the complex masses of the unstable top quark and Higgs boson 
should be taken everywhere in both tree-level and one-loop level calculations. Then the gauge 
invariance is conserved and the singularity poles of propagators for real p 2 are avoided. The 
relevant complex masses are defined as /if = mf — imtTt, ^ 2 H = rrvfj — imH^H, where rat and ran 
are the conventional real masses, Tt and Vh represent the corresponding total widths of the top 
quark and Higgs boson, and the poles of propagators are located at fj, 2 and fi 2 H on the complex 
p 2 plane, respectively. Since the unstable particles are involved in the loops for the 0(a s ) QCD 
corrections, we shall meet the calculations of N point integrals with complex masses. 
The LO total cross section for pp — > W + W~b + X can be expressed as 

°lo(pP ->• W+W'b + X) = 

j dx A dx B [G g / A {x A , fJ,f)G b/B (x B , Lif)vLo(gb ->■ W + W~b, x a xbs, Hf,Hr) + (A++ B)] . 

(2.2) 

There G^/p (i = g,b, P = A,B) represent the PDFs of parton i in proton P, Hf and [i r are 
the factorization and renormalization scales separately, and x A and x b describe the momentum 
fractions of parton (gluon or bottom-quark) in protons A and B respectively. 

II.. 2 QCD NLO corrections 
II.. 2.1 General description 

Our jet recombination procedure is based on the jet algorithm in Ref.[12j. We consider the jet 
events with up to two protojets involved in the final states. The candidate jet events are applied 
by the jet recombination procedure with the following criteria: If the two protojets i and j 
satisfy \J At/ 2 + Acp 2 < R = 0.7 (where A77 and are the differences of pseudorapidity and 
the azimuthal angle between the two protojets), they are merged as one jet and this jet event 
is called a one-jet event, whose four-momentum is defined as Pij,u = Pi t u + Pj,u- The transverse 
momentum and rapidity of the merged jet are calculated from its four- momentum by using the 
formulas p ij>T = \Jpfj, x + Pij, y and Vij = \ In (lg±f|f ) • 

The QCD NLO correction to the pp ->■ W + W~b + X process includes the following compo- 
nents: 
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(i) The QCD one-loop virtual corrections to the partonic process gb — > W + W~b. 

(ii) The contribution of the real gluon emission partonic process gb — > W + W~bg. 

(iii) The contribution of the real light-(anti)quark emission partonic process qb — > W + W~bq, 
where q denotes any of the light-(anti)quarks u, d, c, s, u, d, c, s. 

(iv) The corresponding contributions of the PDF counterterms. 

(v) The additional contributions of the PDFs for subtracting the quasicollinear mass singu- 
larity from massive bottom quark. 

In order to make a comparison, we use two schemes for the collection of the finite contribu- 
tions at 0{a 2 s a 2 w ). In the QCD NLO contribution collection scheme-I, we collect only the five 
components mentioned above for the QCD NLO corrections. In the QCD NLO contribution 
collection scheme-II, in addition to the contributions involved in the collection scheme-I, the 
QCD NLO corrections to the pp — >• W + W~b + X process include the following three additional 
contributions: 

(i) The LO contribution of the (anti)bottom-quark radiative processes qq(gg) — > W + W~bb 
(q = u, d, c, s, b) in the "merged one-jet phase space" (One half contribution belongs to the 
QCD NLO correction of the pp — > W + W~b + X process, another half is to the correction 
of the pp — > W + W~b + X process). 

(ii) The LO contributions of the pp — > bb — > W + W~bb and pp — > bb — >■ W + W~bb processes 
in the "merged one-jet phase space" (The contribution from the former process belongs 
to the pp — > W + W~b + X process, the latter part belongs to the pp — > W + W~b + X 
process). 

(iii) The additional contributions of the PDFs and fragmentation functions (FFs) for subtract- 
ing the quasicollinear mass singularities in the pp — > W + W~bb(bb,bb) processes. 

The so-called "merged one-jet phase space" is defined where the two final b(b) jets are merged 
as one jet in the tagging region of the detector, i.e., \J Ar] 2 + A0 2 < 0.7, and the merged jet 



6 



passes the constraints of pr, j > 25 GeV and \yj\ < 2.5. That is to say the two (anti)bottom 
quarks cannot form two isolated hard jets, and therefore we call the event a one-jet event. We 
notice that all the contributions involved in the QCD NLO contribution collection scheme-II are 

at 0(a 2 s a 2 ew ). 

The dimensional regularization method in D = 4 — 2e dimensions is used to isolate the UV 
and IR singularities. We split each collinear counterterm of the PDF, 5Gi/p(x, /if) (P = proton, 
i = g,u, u,d, d,c, c,s, s), into two parts: the collinear gluon emission part 6G^p° n ^ (x,fif) and the 
collinear light-quark emission part SG^ 9 , p ark \x, fif). The analytical expressions are presented as 
follows. 



5G n 



q{g)/P {x,fi f ) = 5G^° n p \x,fi f ) + 5G^ ( u g yp'(x,fi f ), (q = u,u,d,d,c,c,s,s), (2.3) 



i(quark) 



where 



5G 



(gluon) 

q(g)/P 



(x,n f ) = - 



a s T(l - e) ( 47T/U: 



5G 



(quark) 
q/P 



0,/if) = - 

e 
1 



8G^ k \x, N ) 



2yrr(l -2e) I n) 

a s r(l-e) / 4tt^ 
2vrr(l - 2e) 1 /ij 

a s r(l-e) / 4tt^ 
2tt r(l - 2c) 1 nj 



1 dz 

— p qq(gg)( z ) G Q(g)/p( x / z ^f)' 
dz 

—P qg (z)G g/P (x/z,(j,f), 
1 dz 



z 

d,d,c,c,s,s 



E 



J q=u,u 



P gq (z)G q/P (x/z,ti f ),(2A) 



and the explicit expressions for the splitting functions Pij(z), (ij = qq, qg, gq, gg) can be found 
in Ref.l 



II. .2. 2 Virtual and real emission corrections to gb — > W + W b 

The amplitude at the one-loop level for the partonic process gb — > W + W~b in the SM contains 
the contributions of the self-energy, vertex, box, pentagon and counterterm graphs. In Figj2]the 
four pentagon Feynman diagrams are shown as representatives. 

In order to remove the UV divergences in virtual corrections, we need to renormalize the 
strong coupling constant, the masses and wave functions of the relevant colored particles. In 
our calculations we introduce the following renormalization constants: 



1 + \ 5Z b(t)J ^b(t)i m° b =m b + 5m b , 
IH + SfH, G° = (1 + -5Z g )G li , g° s =g s + 5g s , 



(2.5) 




(1) " (2) " (3) " (4) b 

Figure 2: The pentagon Fynman diagrams for the partonic process gb — > W + W~b. 



where g s denotes the strong coupling constant, m& and [it are bottom-quark mass and top-quark 
complex mass, respectively, and ip^' , ij) t ' R : and denote the fields of bottom quark, top quark 
and gluon, separately. The masses and wave functions of the colored fields are renormalized in 
the on-shell scheme, and the relevant counterterms are expressed as 

,2 



5Z^ R = 


a s {n r ) 
3vr 


bZ^ R = 
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(2.6) 



where Auv = l/o/V — 7_b + ln(47r) and A 77? = 1/em — + ln(47r). 

For the renormalization of the strong coupling constant g s , we adopt the MS scheme at 
the renormalization scale /x r , except that the divergences associated with the massive top- and 
bottom-quark loops are subtracted at zero momentum[14J. Then the counterterm of the strong 
coupling constant g s can be obtained as 



9s 



4tt 



ft 



1 



mt 



Am/ H — In 

2 uv 3 nl 



+ - In 



(2.7) 



3 iv 



|iVy — |, and the numbers of colors and active flavors are taken as N = 3 and 



where (3q 

Nj = 4 respectively. As shown in Eqs. (|2.6p and (|2.7p . the renormalizaion constants SZ g , 5Z^' R , 
dmi) and 5g s /g s contain the terms of a s In ("ir)) which is divergent when m^ — > 0. Therefore, 
the amplitude of the counterterms for the subprocess gb — > W + W~b also contains the quasi- 
collinear mass-singular term. By adding the virtual and real corrections to the contributions of 
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counterterms, the 0{a 2 s a 2 w ) NLO QCD corrected partonic cross section still involves the quasi- 
collinear mass-singular term, which could violate the convergence of the perturbative series. But 
this quasicollinear mass-singular term can be canceled exactly by the additional contributions 
of the PDFs and the FFs for subtracting the quasicollinear mass singularity due to the massive 
bottom quark, as explained in the following subsection. 

Because we use the CMS to deal with the possible top-quark and Higgs-boson resonances, the 
normal one-loop integrals must be continued onto the complex plane. The formulas for calcu- 
lating the IR-divergent integrals with complex internal masses in the dimensional regularization 
scheme are obtained by analytically continuing the expressions in Ref.[15] onto the complex 
plane. The numerical evaluations of IR-safe N point (n = 1,2,3,4,5) integrals with complex 
masses, are implemented by using the expressions analytically continued onto the complex plane 
from those presented in Refs.[16l [XT] [18]. In this way, we created our in-house subroutines to 
isolate analytically the IR singularities in integrals and calculate numerically one-loop integrals 
with complex masses based on the LoopTools-2.4 package[10J[19j. 

The corrections due to the real gluon and (anti) light-quark emission partonic processes, 
KPi) 9[q,q](P2) -> W + (p 3 ) W~(p4,) b(p 5 ) g[q,q](pe)(q = u,d,c,s), are dealt with by using the 
two cutoff phase-space slicing method [15]. The soft IR singularity of the real gluon emission 
subprocess can be isolated by introducing an arbitrary soft cutoff 5 S , to separate the 2 — > 4 
phase space into two regions, Eq < 5 s V§/2 (soft gluon region) and E§ > 5 s V§/2 (hard gluon 
region). The cutoff 5 C decomposes the real hard gluon/ (anti) light-quark emission phase-space 
region into hard collinear {HC) region, S2Q < S c s (where §ij = (pi + Pj) 2 ), and noncollinear 
(HC) region, ,§26 > S c s, in order to isolate the collinear singularity from the IR-safe region. The 
integration over the HC region of phase space can be performed in the four-dimensions by using 
the program based on the VEGAS, a Monte Carlo integrator [20]. Then the cross section for the 
real emission processes can be written as 

Aa R = Acts + Aa H = Aa s + Aa HC + Aaj^. (2.8) 

The UV singularities in the virtual corrections are canceled by the contributions of all the 
related counterterms. Soft and collinear IR singularities are also involved in virtual corrections. 
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After combining the contributions of the real gluon /light-quark emission processes and the PDF 
counterterms SG q r g yp with the virtual contributions together, these IR singularities are exactly 
vanished. These cancelations have been verified numerically in our numerical calculations. 

II.. 2. 3 Subtraction of quasicollinear mass singularity 



Since the bottom-quark mass is retained within the whole calculation, the PDFs of (anti)bottom 
quark and gluon (up to the a s order) contain the large logarithm terms of ^a s In \ • The large 
logarithm terms in the G b / b y P (x, fif) and G g /p(x, yu/) PDFs arise from the (anti)bottom-quark 
emission off the (anti)bottom quark or gluon and gluon emission off the (anti)bottom quark, 
respectively. They are just the QCD NLO counterterms of G b ^ b y P (x, (if) and G g j P {x, (if) which 
are finite with a nonzero bottom-quark mass. The isolation of these quasicollinear mass-singular 
terms are similar to that in the conventional massless parton model approach. 

Analogous to the PDFs of light quarks (q = u,u,d,d,c,c,s,s) in Eq. ([2.3|) . the QCD NLO 
counterterm of G b / b y P (x, fj,f) can be split into the gluon emission part and the b(b) emission 
part. 



SG b(b)/p(. x > M/) = 5G bS ( x ' ^/) + 6G b(l)/p( x > ^/) 



(5(6)) 



6(6) IP 



6(6) /-P 



(2.9) 



The explicit expressions for these two components can be expressed as |214 [22] 



a, f dz 



b(b)/p 



2tt J x z 



P bb(bb)(z) Nn^-21n(l-s)-l 



5G 



(6(6)) 
6(6)/P 



(x,fi f ) = ~ [ T Pb 9(h)^ ( ln ^ ) G g/P (x/z,fi f ), 



int 



G b(b)/p( x / z ^f) 



(2.10) 



where the [...]+ prescription is defined by 



dz[f(z)] + g(z) 



dzf(z)g(z)- f dzf(z)g(l) . 
Jo 



Obviously these two counterterms are proportional to \yj> s In ^ J and finite with nonzero 
bottom-quark mass. 

The QCD NLO counterterm of G g / P (x, fif) is also divided into two components. Among 
them one is the divergent, whose analytical expression is proportional to ^ f^j^fz^ry f ~j^ r ' J 



(2.11) 
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and has been given in Eqs. (|2,3p and (|2.4p . Another component is the additional part, dG^p (x,fif) 
which is proportional to ( a s In — 4 ] , i.e., 



SG^fix^f) = 8Gf /p (x, H ) + 8G® p {x,n f ) 



2vr 



dz ( fJ- } \ 

—Pgb(z) In -2 G b/P (x/z,n f ) 



mt 



(2.12) 



In order to include the contributions of these quasicollinear mass-singular counterterms of 
the PDFs in the total NLO cross section, we replace the PDFs in the expressions for ulo{pP ~~ ^ 
W + W~b + X) and a L o{pp -»• W + W~b + X) (see Eq.Q for a LO (pp -> W + W~b + X) and 
similar one for ctlo(pP —> W + W~b + X)) as 



G b(b)/p( x ^f) — ► G 6(5)/pfo/V)+ 5G 6(S)/£ \ x ,Vf) + SG i { ty P ( x ,Vf), 



G g /p(x,fif) — ► G g /p(x,pf) + SG^Jp 1 ™' (x, nf) + 5G 



i(quark) 

'g/P 1 



x 



(2.13) 



Then we obtain the following five additional contributions of the PDFs: 



(i) 



m 



(iv) 



<J PDF (W + W-bg) 



dxAdxB 



G g/A {x A , H )8G^\x B , H )a L o{,gb -> W+W~b) + (A++ B) 



<j PDF (W + W-bg) 



dx^dxB 



Gg/AixA^f^G^ixB^f^Loigb -> W^Wfe) 



<rjw(W + W~W>) 



dx A dxB 



G b/A (xA,Vf)6G g b ) B (x B ^f)aLo(gb -> + (A « B) 



£Jpbf(^ + VF-66) 



dxAdxB 



G- b/A (xA,Hf)5Gf} B (x B ,lif)&Lo(gb + (4 o B) 

<r P r>j?(W~ , "W~&&) 



dxAdxB 



G g/A (x A ,Hf)SG { b b J B (xB,fif)a L o(gb -> PF+VF 6) 
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+ 



+ G g/A (x A ,ii f )6G&} B (xB,l*f)&Lo(gb -> W + W-~b) + {A^B) 
j dx A dx B G b / A (x A ,fi f )5G^ B (xB,fif)vLo(gb^W + W~b) 

+ G- b/A (x A ,n f )8Gf) B (x B ,fi f )a LO (gb -> PF+^"&) + 5)] . (2.14) 

From the above discussion we can see that the large logarithm in the cross section a{pp — >• 
gg ->■ VF + VF~66+X) is canceled exactly by that in the additional PDF contribution 
shown in Eq. (12.141) . To cancel the quasicollinear mass singularity of thepp — > bb — > W + W~bb+X 
process, both the additional contributions of the PDFs and the FFs should be included, since 
the quasicollinear mass singularity arises not only from the (anti)bottom-quark emission off the 
initial (anti)bottom quark, but also the gluon splitting to the bb pair in the final states. For 
the pp —> qq — >■ W + W~bb + X (q = u,d, c, s) processes there also exists a large logarithm cor- 
responding to the final gluon splitting to the bb pair, which can be absorbed by the additional 
FF contribution. The FF additional contribution to the cross section can be expressed as[22j 



a^DiW+W-bb) 



= -J dx A dx B dz G q / A {x A ,Hf)G q / B (x B ,fi f )a LO (qq W + W g)^P hg (z) [hi ^ 

+(A*+B) , (q = u,d,c,s,b). (2.15) 

II. .3 Cancellation of quasicollinear mass-singularity 

The quasicollinear mass singularity of the real gluon emission process pp —> gb — > W + W~bg + X 
arises from the gluon emission off either the initial or the final bottom quark. The large logarithm 
term corresponding to the gluon emission off the initial bottom quark is canceled by that of 
ap£>F(W + W~bg), while the large logarithm term corresponding to the gluon emission off the 
final bottom quark is canceled by that of the virtual corrections to the pp — > W + W~b + X 
process exactly. For the real light-quark emission processes pp —> qb — )• W + W~bq + X (q = 
u,u,d,d,c,c,s,s), the final bottom quark is noncollinear to the initial bottom quark since a 
bottom jet should be detected in the final states. Therefore, the real light-quark emission 
processes pp — >■ qb — )• W + W~bq + X (q = u,u, d, d, c, c, s, s) do not contain quasicollinear mass 
singularity. 
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m 2 b 



After the cancelation of quasicollinear mass singularity explained above, we get the total 
QCD NLO correction to the pp — > W + W~b+X process by adopting the QCD NLO contribution 
collection scheme-I (see next section), which is free of the contribution from the large logarithm 
term ^a s In £0 , expressed as 

Aa NLO {pp -> W + W~b + X) 

U } d } C,S 

Ao R {pp^ qb^W + W~bq + X)+[A(j v {pp^ gb^W + W~b + X) 

q=u,d,c,s 

+ Aa R (pp^ gb -> W + W~bg + X) + a PDF (W + W-bg)] , (2.16) 

We find that A<jnlo{pP — > W + W~b + X) is convergent when m& — > 0. Analogously, the total 
QCD NLO correction to the pp — > W + W~b + X process, which is the summation of the virtual, 
real corrections and apDp(W + W~bg), is also finite when m& — > 0. 

There are 28 Feynman diagrams for the pp — > bb — > W + W~bb + X process at the 0{a 2 s a1 w ). 
We can find that one (and only one) of the two final bottom quarks is emitted off the initial 
bottom quarks for each of these Feynman diagrams, therefore, the cross section for thepp — > bb — > 
W + W~bb + X process contains the large logarithm term of a s In {j^j ■ This quasicollinear mass 
singular term is canceled exactly by that of <jpp, F (W + W~bb). Similarly, the large logarithm 
terms in a(pp — > bb — > W + W~bb + X) and app, F {W + W~bb) are canceled by each other. Then 
the total cross sections for the pp — > W + W~bb + X and pp — > W + W~bb + X processes defined 
as 

<?tot(pp^ W + W~bb + X) = <j{pp^bb^W + W-bb + X)+a PDF (W + W-bb) 
(J to t(pp^W + W~bb + X) = a(pp^bb^ W + W-bb + X)+<j PDF {W + W~bb), (2.17) 

are convergent when — > 0. 

The quasicollinear mass singularity sources in the pp — > W + W~bb + X process are listed as 
follows: 

(i) pp — > qq — >■ W + W~bb + X (q = u,d, c, s): the internal gluon splitting into 66-pair. 

(ii) pp — >■ gg — >■ W + W~bb + X: the (anti)bottom quark emitting off the initial gluon. 
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(iii) pp — > bb — > W + W bb + X: the (anti)bottom-quark emitting off the initial (anti)bottom 
quark and the internal gluon splitting into bb pair. 

The LO total cross section for the pp — > W + W~bb + X process defined as 

atotipp^W+W-bb + X) = 

[ a ( p P W + W-bb + X)+ a { fJ(W + W-bb) 

q=u,d,c,s 

+ [<r(pp ^bb^ W + W-bb + X) + a ( p^ F (W + W- bb) + af^W^Wbl) 

+ U{PP ^99^ W + W-bb + X)+ a^ F (W + W-bb)] , (2.18) 

which is finite when nib and would be involved in the QCD NLO corrections to the pp — > 
W + W~b + X and pp — > W + W~b + X processes. 

III. Numerical results and discussions 
III..1 Input parameters 

In this work we take one-loop and two-loop running a s in the LO and NLO calculations, 
respectively [23 j . For simplicity we set the factorization scale and the renormalization scale 
being equal (i.e., fi = fij = fi r ) and take n. = fi = rrib/2 + mw in default unless other- 
wise stated. Throughout this paper, we take a eiu (m^) _1 | j^s = 127.925, myy = 80.398 GeV, 
mz = 91.1876 GeV, sin 2 6 W = 1 — {j^^j = 0.222646, and set quark masses as m u = m<i = 
m c = m s = 0, rrif) = 4.2 GeV and nit = 171.2 Gey[23 . The colliding energy in the proton- 
proton center-of-mass system is taken as ^/s = 14 TeV for the LHC. The CKM matrix elements 
are set as 

/ V ud Vus V ub \ I 0.97418 0.22577 \ 
Vckm = V cd V cs V cb = -0.22577 0.97418 . (3.1) 
\ Vtd v ts V th j \ 1 J 

We adopt the CTEQ6L1 and CTEQ6M PDFs in the LO and NLO calculations, respectively. 

The QCD parameter = 215 MeV for the CTEQ6L1 at the LO, and Af^ = 326 MeV for 

the CTEQ6M at the NLO [21]. 

Since we take V t b ~ 1, the decay of the top quark is dominated by the t — > W + b decay 

mode, and the total decay width of the top quark is approximately equal to the decay width of 
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t —> W + b. Neglecting the terms of order m 2 /m 2 , a 2 and (a s /ir)m^ /mf , the width predicted 
in the SM is [25] 



t 



ot ew m\ _ mw\ 2 L + 2m% r \ |\ 2a., f 2n 2 5 



1 — ' 



(3.2) 



3vr \ 3 2 

By taking a ew = ^(ml)^ = 1/127.925 and a s (m 2 ) = 0.1024, we obtain T t = 1.3692 GeF. 
The reasonable physical decay width of the Higgs-boson is obtained by employing the program 
Hdecav |26j . where the partial decay width T(H® — y qq) is calculated up to the (D(a^a ew ). Then 
we obtain T H = 0.2965 x 10~ 2 , 1.704 x 10~ 2 and 0.6511 GeV for m H = 120, 150 and 180 GeV, 
respectively. 

The verification of the independence of the total QCD NLO correction to the pp — > W + W~b+ 
X process on the two cutoffs, 5 S and 6 C , is made. It shows that the total QCD NLO correction 
AcrjvxO) which is the summation of the three-body and four-body cross sections, is independent 
of the two cutoffs within the statistical errors. That independence of the full QCD NLO correc- 
tion to the pp —?• + X process on the cutoffs 5 S and 6 C provides an indirect check for 
the correctness of the calculations. In the further numerical calculations, we fix 5 S = 1 x 10~ 3 
and 5 C = 5 s /50. 

III. .2 Event selection criteria 

In this subsection we present the description of the event selection criteria. After applying the 
jet recombination procedure to the protojet events of the pp — > W + W~b(b) + X up to the QCD 
NLO, we obtain the one-jet events and two-jet events. The one-jet event contains only a b(b) jet 
or a merged jet involving at least one 6(6) quark, while the two-jet event contains two isolated 
hard jets which are 6(6) jet and gluon/light-quark jet, respectively. In the following we present 
detailed criteria for how to treat the one-jet and two-jet events. 

(1) For the one-jet events, we collect the events with the constraints on the jet as pr j > 
25 GeV and \ Vj \ < 2.5. 

(2) For the two-jet events (i.e., W + W~b(b)j (j = g,q and q = u,u, d, d, c, c, s, s) production), 
we treat the hard gluon/light-quark jet (j jet), which is noncollinear to the 6(6) jet, either 
inclusively or exclusively. The inclusive and exclusive two-jet event selection schemes are declared 
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as follows: 

(i) In the inclusive two-jet event selection scheme both one- and two-jet events are included, 
no further restriction is applied on the j jet except for the b(b) jet. The b(b) jet should 
satisfy Pr,b(b) > ^5 GeV and 1 2/^(5) I < 2.5. 

(ii) In the exclusive two-jet event selection scheme, the one-jet events are accepted while 
the two-jet events are rejected. Therefore, besides the W + W~b(b) production events, we 
accept the W + W~b(b)j production events only when the b(b) jet satisfies PtMI) > 25 GeV 
and 1 2/^(5) I < 2.5, and another gluon/light-quark jet (j jet) passes the constraint of either 
PT,j < 25 GeV or \yj\ > 2.5. That means the second j jet is undetectable or there is no 
second separable jet being observed for the W + W~b(b)j production events. 

The cross sections for the pp — > W + W~bb+X, pp — > W + W~bb+X and pp — > W + W~bb+X 
processes combined with the corresponding additional contributions of the PDFs and the FFs 
shown in Eqs. (|2.14p and (|2.15p . are free of the quasicollinear mass singularity induced by the 
a s ln {jJ' 2 f r /fn^j term. Part of these contributions are regarded as the QCD NLO corrections to 
the pp —> W + W~b(b) + X process in the QCD NLO contribution collection scheme-II, where 
the two b{b) protojets are merged as one jet. In Table [2 we present the numerical results 
of the integrated cross sections for the pp — > W + W~bb + X, W + W~bb + X, W + W~bb + X 
processes over three different regions of final state phase space with thh = 120, 150 and 180 GeV, 
respectively. 

III.. 3 Dependence on energy scale 

In Figj3] we present the integrated LO and the QCD NLO corrected cross sections for the 
pp W + W~b{b)+X processes at the LHC, as the functions of the renormalization/factorization 
scale. The full curve is for the LO cross section. The dashed and dotted curves are for the QCD 
NLO corrected cross sections by taking the QCD NLO contribution collection scheme-I, and 
adopting the inclusive and exclusive two-jet event selection schemes, separately. The dash-dot- 
dotted and dash-dotted curves are for the results by taking the QCD NLO contribution collection 
scheme-II, and adopting the inclusive and exclusive two-jet event selection schemes, respectively. 
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We can see from FigJ3]that the curve for the LO cross section is relatively smooth, while the 
QCD NLO corrections in the QCD NLO contribution collection scheme-I make the uncertainty 
on fj, rather large. The QCD NLO corrections in the QCD NLO contribution collection scheme-I 
by using the inclusive ("exclusive") two-jet event selection scheme change the LO cross section by 
a factor of between 1.30 (0.51) and 1.99 (1.60) when [i varies from 0.1/io to 3^o- I n other words, 
the scale uncertainties of the QCD NLO corrected cross sections in the QCD NLO contribution 
collection scheme-I are even worse than that of the LO cross section. This enhancement of 
the uncertainty of the QCD NLO theoretical predictions is partially due to using the 5FS in 
calculations, and in this scheme the perturbative series is not ordered strictly by powers of the a s . 
Figure shows also that the QCD NLO corrections in the QCD NLO contribution collection 
scheme-II by adopting the inclusive ("exclusive") two-jet event selection scheme increase the 
LO cross section by a factor varying from 2.12 (1.59) to 3.54 (2.68) in the plotted range of \i. 
That shows the QCD corrections in the QCD NLO contribution collection scheme-II violate the 
convergence of the perturbative series in some ranges of particularly the range of ^ < 0.5/Uo- 
But we find that the total cross section for the W + W~b{b) production up to the QCD NLO in the 
QCD NLO contribution collection scheme-II demonstrates a weaker scale uncertainty compared 
with the LO one. That means the energy scale uncertainty of the LO cross section can be 
improved by taking the QCD NLO contribution collection scheme-II and adopting either the 
inclusive or the exclusive two-jet event selection schemes. It just reflects that a stable prediction 
of the integrated cross section for the pp — > W + W~b(b) + X process requires not only a veto on 
a second isolated hard jet, but also the inclusion of the QCD NLO corrections contributed by 
the W + W~bb(bb, bb) production events with the two final b(b) quarks being merged as one jet. 
In the further numerical calculations, we take fi = [1q in default of other statements. 

III.. 4 Dependence on Higgs mass 

As indicated in the introduction, the W + W~b{b) production may serve as one of the important 
backgrounds of Higgs-boson production associated with a b(b) jet, if the Higgs-boson mass is 
larger than 2mw- In this subsection we study the behavior of the QCD NLO corrected cross 
section for the pp — > W + W~b + X process as a function of the Higgs mass. According to the 
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Figure 3: The dependence of the LO and NLO QCD corrected total cross sections for the 
pp — > W + W~b + X and pp — > W + W~b + X processes on /x/a*o at the LHC. 
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m H (GeV) 


olo (pb) 


a NLO (inc.)(pb) 


K factor 


a NLO (exc.)(pb) 


K factor 


120 


61.86(2) 


102.84(4) 


1.66 


74.56(2) 


1.21 


150 


61.84(2) 


102.82(4) 


1.66 


74.56(2) 


1.21 


180 


61.96(2) 


102.58(4) 


1.67 


75.18(2) 


1.21 



Table 1: The integrated LO and QCD NLO corrected cross sections and the corresponding 
K factors for the pp —> W + W~b + X and pp — > W + W~b + X processes at the LHC by 
taking the QCD NLO contribution collection scheme-I. 

recent Fermilab report [27J, the SM Higgs-boson mass has been narrowed down to 114 ~ 158 GeV 
and 175 ~ 185 GeV. In our numerical calculations we choose Higgs mass values and their 
corresponding widths being T H (m H = 120 GeV) = 0.2965 x 10" 2 GeV, T H {m H = 150 GeV) = 
1.704 x 10~ 2 GeV and T H (m H = 180 GeV) = 0.6511 GeV separately, as provided in Sec lTTOl 
In Table [1] we present the numerical results of the LO and QCD NLO corrected cross sections, 
and the corresponding K factors by taking the QCD NLO contribution collection scheme-I with 
the Higgs-boson mass being 120, 150 and 180 GeV, respectively. The notations of (inc.) and 
(exc.) refer to the results by adopting the inclusive and exclusive two-jet event selection schemes, 
separately. From this table we find that the LO and QCD NLO corrected cross sections with 
both the two two-jet event selection schemes are insensitive to the Higgs-boson mass when mn 
varies from 120 to 180 GeV . That is because the contributions of the Feynman diagrams, which 
correspond to the H°b production with H° subsequently decaying into W pair, are very small. 

III. .5 pT,b and Pt,w ± distributions 

The LO and QCD NLO corrected distributions of the transverse momentum of the b jet (pr,b) 
and the corresponding K-factors for the pp — > W + W~b + X process at the LHC by taking the 
QCD NLO contribution collection scheme-I, are demonstrated in FigJU The full, dashed and 
dotted curves are for the LO, QCD NLO corrected px,b distributions by adopting the inclusive 
and exclusive two-jet event selection schemes, separately. The corresponding K factors are drawn 
in the attached figure below. These two figures show that the NLO QCD corrections with both 
the inclusive and exclusive two-jet event selection schemes always enhance the LO differential 
cross section daLo/dpT,b, and the NLO QCD corrections to the LO differential cross section 
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Figure 4: The LO, QCD NLO corrected distributions of the transverse momentum of the bottom 
jet (pT,b) an d corresponding K factors for the pp — > W + W~b + X process at the LHC by taking 
fj, = yUo, niH = 120 GeV and adopting the QCD NLO contribution collection scheme-I. The 
corresponding K factors are shown in the figure below. 

with the inclusive scheme are always larger than that with the exclusive scheme. 

We plot the LO and the QCD NLO corrected distributions of the transverse momenta of 
bosons Pt,w+ an d Pt,w- f° r the pp — > W + W~b + X process at the LHC by taking the 
QCD NLO contribution collection scheme-I in Figj5|a) and Figj5]^b), respectively. We use the 
full, dashed and dotted curves describing the LO, QCD NLO corrected distributions with the 
inclusive and exclusive two-jet event selection schemes, separately. We can see from the figures 
that the QCD NLO corrections with both the inclusive and exclusive two-jet event selection 
schemes, always increase the LO differential cross sections daLo/dpT,w+ an d daLO / dpT,w- , 
and the enhancement due to the QCD NLO corrections to daLo/dpT,w with the inclusive two- 
jet event selection scheme is larger than that with the exclusive two-jet event selection scheme. 
From the figures we see also that the distributions of the transverse momenta of W + and W~ 
bosons are different. 
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Figure 5: (a) The LO, QCD NLO corrected distributions of the transverse momentum of W + 
boson and the corresponding K factors for the pp — > W + W~b + X process at the LHC. (b) 
The LO, NLO QCD corrected distributions of the transverse momentum of W~ boson and the 
corresponding K factors. In these two figures we take \i = /xo, m// = 120 GeV and adopt the 
QCD NLO contribution collection scheme-I. 

III.. 6 Invariant mass mryv+b) distribution 

We plot the distributions of the invariant mass of the (W + b) pair, mn^+f^, at the LO and the 
QCD NLO in Figj6]with fi = /iq and run = 120 GeV, where the NLO distributions of m^+^ are 
presented by taking the QCD NLO contribution collection scheme-I and adopting the inclusive 
and exclusive two-jet event selection schemes separately. We can see from the figure that the 
LO and QCD NLO corrected differential cross sections, daLo/dfn(w+b) an d daNLo/dm^ w + b ^, 
have peaks in the vicinity of uiryy+fy = m%. It demonstrates that there is a large contribution 
part for the process pp —> W + W~b + X , which comes from the single top production associated 
with a W~ boson (pp — )■ bg — > W~t) followed by the subsequential decay of the top quark to 
W + b. That is to say the main contribution to the pp —> W + W~b + X process originates from 
the Feynman diagrams of Figs{T](4,6) & t the LO, and their related one-loop diagrams at the 
NLO. Here we can see that in the invariant mass range around rnryy+b) ~ 171 GeV the QCD 
NLO corrections with both the inclusive and exclusive two-jet event selection schemes strongly 
enhance the LO invariant mass distribution da lo I 'dmrw+q, and the QCD NLO correction with 
the inclusive scheme enhances the LO distribution more heavily than that with the exclusive 
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Figure 6: The distributions of the invariant mass of the bottom jet and W + boson, mny+fy, for 
the pp — > W + W~b + X process at the LHC, where we take [i = no, mn = 120 GeV and adopt 
the QCD NLO contribution collection scheme-I. 

scheme, particularly around the resonance peak. 

III. .7 Cross sections for the pp — > W + W~bb, W + W~bb, W + W~bb processes 

In Table [2] we list the numerical results of the integrated cross sections for the three pro- 
cesses pp —7- W + W~bb and pp — > W + W~bb, bb with different Higgs-boson mass values. These 
three W + W~bb, W + W~bb and W + W~bb production processes are induced by the gg(qq) — > 
W + W~bb, bb — > W + W~bb and bb — > W + W~bb partonic processes separately. The data in 
columns (I), (II) and (III) are the integrated cross sections over three different phase-space 
regions, respectively. 

In phase-space region (I), the two b(b) jets are isolated (i.e., the two b(b) jets satisfy sj Arj 2 + Acfi 2 > 
0.7), and one b(b) jet is detectable satisfying the constraints oip Tb ^ > 25 GeV and {y^l < 2.5 
while the other 6(6) jet is undetectable, i.e., satisfying PT,b(l) < 25 GeV or 1 2/£,(fo) I > 2.5. 

In phase-space region (II), the two 6(6) jets are merged as one jet (i.e., \j Arj 2 + Acj) 2 < 0.7), 
and the merged jet passes the constraints of ptj > 25 GeV and \yj\ < 2.5. That is just the 
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m H (GeY) 


WWbb(bb,bb){l)(pb) 


WWbb(bb, 66)(II)(pb) 


WWbb(bb, 66)(III)(pb) 


120 


141.66(6) 


27.15(2) 


603.2(1) 


150 


141.54(6) 


27.15(2) 


603.5(1) 


180 


141.28(6) 


27.15(2) 


603.3(1) 



Table 2: The integrated cross sections for the three processes pp —> 
W + W~bb, W + W~bb, W + W~bb over three different phase-space regions (I), (II) 
and (III) with fj, = fiQ, ran = 120, 150 and 180 GeV, respectively. 

additional QCD NLO contribution when the QCD NLO contribution collection scheme-II is 
adopted. 

In phase-space region (III), the two b(b) jets are separately detectable. This means that 
a/ Arj 2 + Acf) 2 > 0.7 and the two b(b) jets separately satisfy the constraints of PtmI) > 25 GeV 
and 12/6(5)1 < 2 - 5 - 

Analogy of the data in Table [H Table [2] shows that the integrated cross sections for the 
processes pp — > W + W~~bb(bb,bb) over phase-space regions (I), (II) and (III) are insensitive 
to the Higgs-boson mass. That means the contributions from the Feynman diagrams, which 
correspond to the H°b production followed by the decay of H° — > W + W~, are rather small. 
We also find from Table [2] that the integrated cross sections for the pp —> W + W~bb, pp —> 
W + W~bb and pp —> W + W~bb processes over the three phase-space regions are all rather large, 
particularly phase-space region (III). Our calculation demonstrates that the contribution from 
the pp —> W + W~bb + X process is the largest one, which has a 2 ~ 3 order larger contribution to 
the cross section than the other two processes. The large enhancement of this process with the 
final bb pair is due to the large gluon luminosity in PDF at the LHC and the overcompensation 
of the top quark resonances in the pp — )• tt — > W + W~bb + X process. 

IV. Summary 

In this paper we calculate the full QCD NLO corrections to the W-pair production in associa- 
tion with a massive (anti)bottom jet in the SM at the LHC, which is an important background 
of the single top production and can be used to search for new physics beyond the SM by 
using the five-flavor scheme. The dependence of the integrated cross sections on the factoriza- 
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tion/renormalization scale and the Higgs mass is studied. We also investigate the QCD NLO 
corrections to the distributions of the transverse momenta (7^7) of the final particles. Our nu- 
merical results show that by taking the QCD NLO contribution collection scheme-I and adopting 
the inclusive (exclusive) two-jet event selection scheme, the K factor of the total cross section 
for the pp — > W + W~b(b) + X process at the LHC can reach 1.66 (1.21) with fj, = /jLq. We find 
that the LO integrated and differential cross sections are modified by the QCD NLO radiative 
corrections obviously, and the QCD NLO corrections to the W + W~b(b) production process with 
the QCD NLO contribution collection scheme-I can not reduce the scale uncertainty of the LO 
cross section for the pp — > W + W~b(b) + X process at the LHC by adopting either the inclusive 
or exclusive two-jet event selection schemes. The stable prediction for the integrated cross sec- 
tion requires not only a veto on a second isolated hard jet, but also the inclusion of the QCD 
NLO correction contributed by the W + W~ bb(bb, bb) production events with the two final b(b) 
quarks being merged as one jet. 
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